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CSALbstract 

cyp discuss magnetism in spinor quantum gases theoretically 
^ild experimentally with emphasis on temporal dynamics of 
^he spinor order parameter in the presence of an external 
Gnagnetic field. In a simple coupled Gross-Pitaevskii picture 

■Mve observe a dramatic suppression of spin dynamics due to 

4_yuadratic Zeeman " dephasing" . In view of an inhomogeneous 
density profile of the trapped condensate we present evidence 

"t^f spatial variations of spin dynamics. In addition we study 

rtpinor quantum gases as a model system for thermodynamics 

|t)f Bose-Einstein condensation. As a particular example we 

■"C^Eesent measurements on condensate magnetisation due to 
Qhe interaction with a thermal bath. 

o 

1 Introduction 

.'Zn^he field of cold quantum gases has seen a rapid growth 

Quince the first realisation of Bose-Einstein condensation in 

T^ijute atomic gases in 1995 Ql |2 El accompanied with the 

^^^^evelopment of a broad range of tools for the detailed con- 

^^r'ol of these systems. Single component Bose-Einstein con- 

^ildensates have evolved into a fundamental model system 

'"^NSll.owing many intriguing phenomena. For an overview see 

t^.g. HElElEllHlElITniin]. In contrast to these systems 

C^ith a scalar order parameter spinor Bose-Einstein conden- 

I sates offer spin as a new degree of freedom and are conse- 

' ^uently represented by a vector order parameter. In addition 

^o being mixtures of different bosonic species the different 

fCqmponents in these multicomponcnt quantum systems are 

^oupled and can exchange particles. This makes spinor Bosc- 

^jilinstein condensates unique systems, which on the one hand 

ossess intrinsic magnetic properties and on the other hand 

ive access to well controlled Bose-Einstein thermodynamics 

with adjustable heat and particle bath. 

Magnetism in degenerate quantum gases offers new regimes 
for studies of collective spin phenomena [T^[T^[T^[T^[T?)l[r7l 
[THlIiniEniinil23linil21i3l2Sll23llHland opens new per- 
spectives in view of the closely related entangled spin systems 
in atomic quantum gases, which show intriguing prospects for 
quantum optics and quantum computation ^29 30 31 321l38j . 
So far studies concentrated on the magnetic properties 
of spin 1 ultracold quantum gases in optically trapped 
23Na 122 ESI El ESI and recently in 87Rb ^1 1^ EHl EHl, 
where also the intrinsically more complex F=2 spin state be- 
came accessible ^7\ ES] . 
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Systems closely related to these spinor condensates are ef- 
fective spin- 1/2 systems realized by radiofrequency coupling 
of two hyperfine states in ®^Rb ^J EDI l^i in which spin- 
waves I^S] and decoherence effects were observed [5^. 

In this paper we will concentrate on F=l and F=2 spinor 
condensates of ^^Rb in two limits. First we investigate the co- 
herent spinor evolution of a trapped ensemble at zero temper- 
ature in the presence of a homogeneous magnetic field, where 
we find suppression of spin dynamics due to the quadratic 
Zeeman effect as well as a spatial depedence of the dynamics. 
The other limit is thermally dominated spin dynamics at tem- 
peratures close to Tc, where a significant fraction of the atoms 
occupies the normal component. Spinor gases in this regime 
can act as a versatile model system for thermodynamics with 
tunable heat and particle bath as recently demonstrated with 
a constant temperature Bose-Einstein phase transition |37j . 
In this paper we will present new data on thermally induced 
condensate magnetisation as another intriguing example of 
spinor thermodynamics. 

2 Spinor condensates at T=0 

The theory presented in this section is based on a mean field 
approach, in extension of the very successful treatment of 
single component Bose-Einstein condensates. The basic two- 
particle interactions are represented by a density and spin- 
composition dependent average energy shift. This approach 
has first been developed for F=l systems |12[ [T^ and was 
later extended to F=2 systems ^| E] . 

For typical experimental parameters the mean field shifts 
connected to collisions in different spin channels dominate 
magnetic dipole dipolc interactions by at least one order of 
magnitude. In the following analysis magnetic dipole dipole 
interactions will thus be neglected. The intrinsic dynamics of 
a spinor condensate is determined by a pairwise interaction 
potential pHIT^ : 
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Here a/ denotes the s-wave scattering length for a collision 
channel of two particles whose single spins F are combined to 
give the total spin /, Pf is the corresponding projection oper- 
ator onto total spin / and m is the mass of a single atom. Due 
to Bose symmetry only even total spin channels (e.g. ao, 02, 04 
for F ~ 2) are involved with the maximum total spin given by 
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/ = 2F. Making use of the relation ( -Fi • F2 i — z^ f 

with Xf — ■^ [/(/ + 1) — 2F(i^ + 1)] the projection operators 
can be replaced by spin expectation values |12| . 

In the following we will concentrate on the case F = I for 
the theoretical considerations in order to point out some im- 
portant aspects of spin dynamics, which are straightforward 
to extend to the F — 2 case. 

In second quantised form the Hamiltonian for a i^ = 1 
system at zero magnetic field is given by 12 : 
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In this expression V'a(r) is the field annihiliation operator for 
an atom in state mp = a at point f and Vext is the trap- 
ping potential. The spin-independent mean-field interaction 
is parameterised by go = ^^^ x 2B2±Sio^ ^j^g spin-dependent 
mean-field responsible for the systems magnetic properties is 

'-^^^ and the 
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characterised by the parameter ^2 = 
coupling between different states is determined by the spin- 
matrices: 
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From the spin matrices it follows, that for F=l spinor conden- 
sates the only coupling process is between states with mp = 
and mp = ±1 and that the total spin projection is preserved 
by the interaction. It is important to emphasise this again to 
demonstrate that the total spin in a finite atomic quantum 
gas system is conserved, in contrast to a homogeneous infinite 
system and in contrast to many condensed matter systems. 
Therefore the actual ground state of the system depends on 
the initial magnetisation. 

For macroscopically occupied Bose-systems at T = it is 
common to replace the field annihilation operators by their 
expectation value, i.e. ipa{r,t) = ('0a(r, t)), which for spinor 
condensates is conveniently expressed as 52] 



Mr,t) = ^/Mr^e'*^''''\a{r,t). 



(5) 



Here n(r, t) is the condensate density, 4>{f, t) a phase and 
C(r, t) = (C+i, Coj C-i)'^ is a normalized spinor with ^t . ^ = x^ 
Using lOl, neglecting the density dependence on the spin 
state and adding the effect of a weak magnetic field one gets 
the following system of differential equations for the evolution 
of an F=l spinor condensate: 
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Figure 1: Absorption images and vertically summed cross 
sections of a spinor condensate after 10 ms evolution start- 
ing with the initially prepared state F = 2, mp = 0. The 
images were taken after a combined time of flight and Stern 
Gerlach separation, leading to a rapid radial expansion and 
a separation of the different spinor components. The graphs 
clearly show a reduced width in the axial distribution of the 
mp = ±1 components created by spin dynamics in the high 
density center of the mp — component. 
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From these equations it follows, that spin dynamics (repre- 
sented by the terms containing T) is proportional to g2n{r, t), 
i.e. the local density. This has important consequences for 
the usual experimental situation of trapped samples, having 
an inhomogeneous density distribution. Two regimes can be 
identified: 1. the spin dynamics rates are higher than one 
or more trapping frequencies and 2. spin dynamics is slow 
compared to trap dynamics. In the first case there will be 
a significant coupling between spin dynamics and motional 
dynamics, while the second case is comparable to the homo- 
geneous density case. As was shown for the first time in |57] 
^^Rb has the fascinating properties that it offers fast spin 
dynamics in the F=2 state (on the order of a few ms) and 
slow spin dynamics in the F=l state (on the order of a few 
s). Therefore also different spatial regimes can be thought of 
with ^^Rb spinor condensates. 

Fig. n shows the initial evolution of an F=2 *^Rb spinor 



condensate prepared in the mp = spin state with a central 
density n « 4 • 10^'*cm~^. In this case spin dynamics takes 
place with timescales on the order of a few ms |22], faster 
than the axial trap dynamics characterised by a frequency of 
« 21 Hz. The other trapping frequencies of ~ 155 Hz and 
» 890 Hz imply motional dynamics faster than the spin evo- 
lution and are thus neglected in the following. The image in 
Fig. n was taken after switching off the trapping potential and 
31 ms time of flight during which a magnetic field gradient was 
applied for 5 ms to achieve a Stern-Gerlach separation of the 
spinor components. The axial trap direction was horizontal, 
along which only weak mean field induced expansion takes 
place during time of flight. The vertically summed horizontal 
cross sections shown in Fig. ^ thus approximately represent 
the axial distribution of the trapped sample. The images and 
graphs clearly demonstrate the density dependence of spin 
dynamics, as the mp = ±1 spin states show a smaller width 
than the "mother" mp — component. This is due to the 
parabolic density profile of the "mother" component, which 
implies a fast spin state conversion in the center of the trap. 
Very recently a spatial dependence in ^^Rb (F=2) condensate 
spin dynamics was also observed in 35 . Further investiga- 
tions of the coupled spin and spatial dynamics were limited 
in our case by the finite optical resolution and will be subject 
to future investigations. 

Interestingly external magnetic fields are a further addi- 
tional parameter to control spin dynamics. Whereas most 
theoretical work so far concentrated on the physics at B=0, 
experiments |2I HZI I28L OS] clearly demonstrated the impor- 
tance of extermal field influences. Magnetic flelds can com- 
pletely hinder spin dynamics or on the other hand strongly 
stimulate dynamics. In the following we will mainly concen- 
trate on a discussion of spin dynamics suppression due to the 
quadratic Zeeman effect. Very recently the suppression of 
spin dynamics in F=2 spinor condensates was experimentally 
observed in [3S]. The following discussion can explain these 
observations in their main part. 

The magnetic field enters equations iQ via the linear Zee- 
man shift as well as the quadratic Zeeman shift with the 



coupling constants p = gp^sB and q = — 
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representing the hyperfine splitting. As we will see later the 
relatively large linear Zeeman energy {q/B w fc^ • 34/iK/G) 
does not infiuence the spin dynamics, which in fact is due to 
total spin conservation. This is not true for the quadratic 
Zeeman effect! For typical experimental conditions with off- 
set fields of several 100 mG the quadratic Zeeman energy 
{\q\/B'^ « ks ■ 3.5nK/G^) can reach values comparable to 
the intrinsic spin coupling (on the order of ks times one nK 
for typical condensate densities of a few 10^"* cm^'^). 

In order to extract the basic influence of the quadratic 
Zeeman effect on spin dynamics, we assume a homogeneous 
case with constant n{f,t) — n and no spatial spin variation 
C{r,t) — C{t)- With these assumptions equation jT)) reads 
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in^^Co = 52ri(C;iC+iCo + CiCoC-i+2Co*C+iC-i) 

-49C0 , (8) 

*?i^c-i - 32ri(c-iC-iC-i-c;iC+iC-i + CoCoC-i 
+c;iCo)+K-i-3<zC-i- 

Using a simple change of variables 

A±i = C±i exp{Ti{p±3q)t/ti) 

and (9) 

Aq = Co exp{-'iAqt/h) , 

the linear Zeeman dependence is removed from the equations 
and the quadratic Zeeman effect enters in a more symmetric 
way: 
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In these equations derived here spin exchange is described by 
the terms with exponentials, while the other terms represent 
the mean field phase evolution. Fig. [2 shows the relative 
spinor occupations and the relative phase (/)+i — 2(f>Q + 0_i 
as a result of a numerical simulation of equations (|10|l for 
different magnetic offset fields. The spin component phases 
are given by (/>i = arg(Ai) with the phase 0+1 — 2(/>o-f 0_i is the 
relevant phase for the evolution of the mp = component. 
The initial conditions were chosen symmetric with density 
n = 4-10i''cm~3, |Ao| =0.9, |A±i| =0.05 and 0j(i = 0) = 0. 

The numerical simulation clearly shows that the quadratic 
Zeeman effect strongly suppresses spin dynamics if it is larger 
than the spin dependent mean field shifts. For *^Rb in F=l 
and at densities of a few times 10^"* cm^^ this suppression 
becomes relevant at magnetic fields of a few 100 mG. 

The suppression of spin dynamics at high magnetic fields 
can also be directly deduced analysing the expression for Aq, 
which according to equation (|10|l at high magnetic fields will 
approximately evolve due to the rapidly changing exponential 
giving: 
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If g 3> 52'^ there will be nearly no change in the occupation 
of spin states as predicted by the numerical simulation. 

In summary (in the absence of a field gradient) the linear 
Zeeman effect can be neglected for investigations on spin dy- 
namics in spinor Bose condensates due to the always symmet- 
ric exchange of Zeeman energy (Fig.|2^), which is fundamen- 
tally caused by spin conservation. In contrast we found that 
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Figure 2: Numerical simulation of equations H1U|) for different 
magnetic field strengths. The graphs on the left show the 
evolution of the population of spin state mp = (red) and 
of spin states m^? = ±1 (black). The graphs on the right 
show the evolution of the relative phase corresponding to the 
interaction 2|0) ^^ |1) + | — 1). A strong reduction of the 
amplitude of the oscillation in the spin populations for higher 
magnetic fields is clearly visible. 
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Figure 3: Schematic view of F=l (|1) + | - 1) ^ 2|0)) spin 
dynamics with additional magnetic field. The different mp 
states are labeled by and ±1. (a) The total linear Zeeman 
energy is conserved in spin dynamics as the energy gain in 
one component is lost in the other, (b) The quadratic Zeeman 
shifts (shown with the linear Zeeman contribution subtracted 
and rescaled as typically AEq ^ AEl) however lead to an 
energy imbalance in spin dynamics. 



spin dynamics can be significantly altered by the quadratic 
Zeeman effect, where an additive energy exchange (Fig. ^p) 
leads to a "dephasing" of the spin components ultimately 
stopping spin dynamics at high magnetic fields. One could 
say that high external magnetic fields "pin" the spin to its 
value. Intrinsic spin dynamics can be observed up to mag- 
netic fields for which q{B) ~ g2n typically corresponds to 
fields of a few 100 mG. We want to emphasize that the block- 
ing of spin dynamics is solely due to the quadratic Zeeman 
effect and docs not follow energetical considerations. Indeed 
with this blocking effect we can explain the experimentally 
observed high magnetic field suppression of spin dynamics 
even when it is leading to an energetically lower state f35] . 

Furthermore we found evidence of spatially varying spin 
dynamics in trapped spinor condensates with inhomogeneous 
density. This effect will lead to complex coupled dynamics 
of spatial and spin degrees of freedom to be investigated in 
future experiments. 

3 Thermodynamics with spinor con- 
densates 

Finite temperature effects in Bose-Einstein condensates rep- 
resent an active area of research, which is still relatively un- 
explored due to its theoretical complexity and experimental 
challenges. In theory sophisticated methods have been devel- 
oped to reduce the complexity of simulations such that mod- 
elling complex phenomena seems feasible [HHl EHl QHI ^1 Q2] ■ 
Major tests for these models consisted in the interpretation of 
early experiments on damping of single component conden- 
sate excitations in the presence of a normal component and 
on condensate formation. 

Spinor condensates offer a novel approach to well controlled 
Bose-Einstein thermodynamics. As a first aspect they are 
multicomponent systems such that a thermal bath for one 
component can easily be created by tailoring the other com- 
ponent (s). This aspect is widely used in sympathetic cooling 
experiments with multi-species mixtures. Spinor multicom- 
ponent systems add the essential aspect of particle exchange 
between the components, which is required to complete ther- 
modynamics. The particle exchange takes place due to in- 
trinsic intcrparticle interactions but can also be experimen- 
tally controlled via additional external electromagnetic fields. 
These aspects and the relevant interactions with respect to 
thermodynamics in spinor quantum gases are shown schemat- 
ically in Fig. ^ 

An important point in the study of spinor systems is con- 
nected to the coherence between different spin components. 
Coherent spin mixtures, i.e. mixtures in which each single 
atom (in the normal component as well as in the condensate 
fraction) is in the same quantum superposition of spin states, 
are effectively single component quantum gases, which have 
to be contrasted to incoherent spin mixtures, where different 
spin states represent different species gases. For example if 
a F=l spinor gas in the first case would be described by N 
particles in a spin superposition state cj) = a\l) +/3|0) +7I — 1) 
then the incoherent state would be given by a mixture of three 
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Figure 4: Schematic view of thermodynamics in spinor quan- 
tum gases for the example of a spin F=l system. In this 
example the states with mF = 1 = |1) and mp — —1 = | — 1) 
act as a heat bath and particle reservoir for the state mp — 
= |0). Heat (energy) is exchanged by elastic collision pro- 
cesses involving atoms in the normal component. Particles 
are exchanged in spin-changing collisions via the interaction 
1 1) -I- I — 1) ^^ 2|0). An important aspect of particle ex- 
change lies in the density dependence of spin dynamics, which 
nearly exclusisively takes place in the dense condensate frac- 
tion. Particle exchange thus involves only small energy trans- 
fer and does practically not contribute to thermalization pro- 
cesses. Spinor condensates thus allow to create systems with 
independently tunable particle and heat exchange. 



gases, one with Ni = |ap particles in the |1) state, one with 
^0 = |/3p particles in the |0) state and one with 7V„i = I7P 
particles in the | — 1) state. 

Indeed incoherent spin mixtures are in some cases of high 
experimental importance, e.g. for the conversion of two spin 
state fermion mixtures to a molecular Bose gas with a Fesh- 
bach resonance. This was nicely demonstrated and explained 
in j43| for the preparation of a spin state mixture in ^Li. As 
an another example the distinction between coherent and in- 
coherent spin superpositions is crucial to the understanding 
of the recently demonstrated decoherence driven cooling |35) 
in a quasi spin 1/2 system. 

The evolution of F=l and F=2 spinor condensates dis- 
cussed in this paper can be tuned in between the regimes 
of coherent and incoherent evolution by adapting the pa- 
rameters temperature, density and possibly external radiofre- 
quency coupling. An intriguing example for mostly coherent 
evolution is the observation of spinor oscillations [271 I28L 15^ . 
The incoherent limit was recently reached in a thermaliza- 
tion dominated regime with the demonstration of constant 
temperature Bose-Einstein condensation ^J in F=l spinor 
condensates with significant occupation in the normal com- 
ponent. 

In the following we will further concentrate on the ther- 
malization dominated regime in F=l ^^Rb and investigate 
the evolution of an initially prepared |1) -I- |0) mixture (see 
Fig.EJ. 

Due to total spin conservation the only spin dynamics is the 
coupling 2|0) ^^ |1)-|-| — 1), which in this case initially leads to 
the depletion of the |0) state in favor of the (initially already 
populated) |1) and the (initially empty) | — 1) states. As spin 
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Figure 5: Absorption images and cross sections showing the 
temporal evolution of a spinor condensate prepared in the 
states mp — 1 and mp = with a significant fraction of 
atoms in the normal cloud (the initial population in mp = — 1 
is due to slight preparation imperfections). 



dynamics is mostly occurring in the condensate fraction, low 
energy atoms are added to the |1) and | — 1) states, which in 
the case of the |1) state just add to the condensate fraction 
(the particle number in the normal component of this state is 
saturated for the given temperature). The case of the newly 
populated | — 1) state is however significantly different, as 
this state does not yet possess a normal component. The 
low energy atoms in this state quickly (on a shorter timescale 
- on the order of 50 ms than spin dynamics - on the order 
of seconds) thermalise with the normal component atoms of 
the other spin states. This leads to a slow buildup of the 
I — 1) normal component and at the same time an increase in 
condensate fraction in the |1) component, while the |0) state 
condensate fraction decreases. 

An interesting point is that this process leads to a slight 
decrease in temperature, as it uses energy from the existing 
normal components to thermalise the low energy atoms en- 
tering the I — 1) state from the |0) condensate fraction. This 
temperature decrease at the expense of total condensate frac- 
tion is similar to decoherence driven cooling observed in quasi 
spin 1/2 systems PEl- 

We want to emphasise that under typical experimental con- 
ditions the thermal energy corresponds to roughly ks x 300 uK 
and is thus more than an order of magnitude larger than the 
spin dependent mean-field shifts of roughly ks x 10 nK re- 
sponsible for spin dynamics. This directly implies that in 
thermal equilibrium the normal components of different spin 
states will have equal population (if there are sufficiently 
many atoms available in each spin component). For the case 
discussed in this paper the | — 1) normal component will grow 
until either the |0) condensate fraction is completely depleted 
or it reaches its saturated occupation for the given tempera- 
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Figure 6: Schematic view of condensate magnetisation in 
thermally dominated spin dynamics. The top row shows the 
initially prepared state. Spin dynamics is transferring popu- 
lation from the |0) condensate fraction to the ±1) condensate 
fractions (second row). Due to fast thermalisation and the ab- 
sence of a I — 1) normal component, the new | — 1) atoms start 
to populate the normal component, while the new | + 1) just 
add to the corresponding condensate fraction (third row) . In 
the end there will be thermal equilibrium with equally pop- 
ulated normal components (zero total spin) and a stronger 
than initial condensate fraction magnetisation. 



ture, i.e. the same occupation as the other spin state normal 
components (which are saturated, as there exists a conden- 
sate fraction in these components). Only in the second case 
a condensate fraction will build up in the | — 1) state, which 
now will be determined not by the thermal energy scale but 
by spin dynamics |37| . 

In any case this process will tend towards a total zero spin 
in the normal components (equal occupation) and thus shift 
the total spin of the condensate fractions towards more pos- 
itive values. The strongest magnetisation of the condensate 
fraction occurs if the initial population of the |0) condensate 
fraction does not suffice to saturate the | — 1) normal com- 
ponent occupation via spin dynamics. In this case spin dy- 
namics stops after the |0) condensate fraction is depleted and 
only a I + 1) condensate fraction remains, i.e. the condensate 
fraction is fully magnetised. 

The principle mechanisms are again summarized in Fig. 
The population of the | — 1) condensate part thermalises and 
fills up a | — 1) normal component. Thus the normal conponent 
total spin finally adds up to zero. Due to spin conservation, 
the condensate spin has to increase, which is reflected in a 
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Figure 7: Experimental data versus a simulation for thermally 
dominated spinor dynamics for an initial preparation of the 
sample in mp = 1 and nip = 0. The graph shows the average 
magnetisation of atoms in the condensate fraction, atoms in 
the normal fraction and in total. 



higher | + 1) condensate fraction population. This process is 
clearly reflected in the experimental data presented in Fig. [3 
The data is well reproduced by a numerical simulation based 
on a simple rate equation model, presented in detail in |37|. 
The slight decrease in the average spin for the total ensemble 
is due to the fact that the trap losses are dominated by three 
body collisions, predominantly occurring in the magnetised 
condensate fraction. 

In conclusion in this paper we have presented investiga- 
tions on spatial variations, influence of magnetic fields and 
high temperature as fundamental and new aspects in spinor 
dynamics. We found that fast spin dynamics in inhomo- 
geneous (trapped) ensembles leads to spatial effects which 
promises new complex coupled spatial and spin dynamics. 
We have shown, that the dominant magnetic field influence 
stems from the quadratic Zeeman effect, limiting the offset 
fields up to which spinor dynamics can be observed for typi- 
cal experimental conditions to a few 100 mG. Furthermore we 
investigated the regime of finite temperature spinor dynam- 
ics considering the example of condensate magnetisation in 
favour of an equalised spin distribution in the normal compo- 
nent. This work demonstrates the versatility and complexity 
of spinor Bose-Einstein condensates and paves the way for a 
broad range of future investigations. 

We acknowledge support in SPP1116 of the Deutsche 
Forschungsgemeinschaft . 
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